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THE  COHERENT  FLAME  MODEL  FOR  TURBULENT  CHEMICAL  REACTIONS 


Abstract 


A description  of  the  turbulent  diffusion  flame  is  proposed  in  which 
the  flame  structure  is  composed  of  a distribution  of  laminar  diffusion 
flame  elements,  whose  thickness  is  small  in  comparison  with  the  large 
eddies.  These  elements  retain  their  identity  during  the  flame  development; 
they  are  strained  in  their  own  plane  by  the  gas  motion,  a process  that  not 
only  extends  their  surface  area,  but  also  establishes  the  rate  at  which  a 
flame  element  consumes  the  reactants.  Where  this  flame  stretching  process 
has  produced  a high  flame  surface  density,  the  flame  area  per  unit  volume, 
adjacent  flame  elements  may  consume  the  intervening  reactant,  thereby 
annihilating  both  flame  segments.  This  is  the  flame  shortening  mechanism 
which,  in  balance  with  the  flame  stretching  process,  establishes  the  local 
level  of  the  flame  density.  The  consumption  rate  of  reactant  is  then  given 
simply  by  the  product  of  the  local  flame  density  and  the  reactant  consump- 
tion rate  per  unit  area  of  flame  surface.  The  proposed  description  permits 
a rather  complete  separation  of  the  turbulent  flow  structure,  on  one  hand, 
and  the  flame  structure,  on  the  other,  and  in  this  manner  permits  the 
treatment  of  reactions  with  complex  chemistry  with  a minimum  of  added 
labor.  The  structure  of  the  strained  laminar  diffusion  flame  may  be 
determined  by  analysis,  numerical  computation,  and  by  experiment  without 
significant  change  to  the  model. 

The  flame  density  and  the  mass  fractions  of  reactant  are  described 
by  non-linear  diffusion  equations  in  which  those  equations  for  the  reactants 
each  contain  a consumption  or  production  term  proportional  to  the  local 
flame  density.  The  flame  density  equation  contains  a production  term 
associated  with  flame  surface  stretching  and  a consumption  term  describing 
the  flame  shortening  by  mutual  annihilation.  Each  of  the  equations  contains 
a turbulent  diffusion  term  utilizing  a scalar  diffusivity.  The  model  of 
inhomogeneous  turbulence,  proposed  by  Saffman,  completes  the  description 
of  the  problem  and  couples  with  the  flame  and  composition  equations  to 
determine  the  velocity  distribution  and  the  turbulent  diffusivity.  A 
single  additional  universal  constant,  over  those  appearing  in  Saffman' s 
model,  is  required  in  the  model  equations  for  the  flame. 

The  coherent  flame  model  has  been  applied  to  diffusion  flame  structure 
in  the  mixing  region  between  two  streams  and  predicts  correctly  the  result 
that  the  reactant  consumption  per  unit  length  of  flame  is  independent  of 
the  distance  from  the  initiation  of  mixing.  In  this  example  which  is 
carried  out  for  small  density  changes,  both  the  fluid  mechanical  and  flame 
variables  possess  similarity  solutions. 

The  coherent  flame  model  is  also  applied  to  the  turbulent  fuel  jet 
which  clearly  does  not  have  a similarity  solution  simply  because  the  finite 
mass  flow  of  fuel  is  eventually  consumed.  The  problem  is  solved  utilizing 
an  integral  technique  and  numerical  integration  of  the  resulting  differential 
equations.  The  model  predicts  the  flame  length  and  superficial  comparison 
with  experiments  suggest  a value  for  the  single  universal  constant.  The 
theory  correctly  predicts  the  change  of  flame  length  with  changes  in  stoich- 
iometric ratio  for  the  chemical  reaction. 
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1.  INTRODUCTION 


Some  of  the  earliest  considerations  of  turbulent  combustion  processes, 
e.g.,  Damkohler/1)  Shchelkin/2'1  included  the  suggestion  that  a turbulent 
combustion  field  consists  of  a collection  of  laminar  flame  surfaces  in  which 
the  laminar  flame  structure  retains  its  identity  while  being  distorted  by  the 
turbulent  motion.  This  conjecture  found  additional  support  a decade  later  in 
the  work  of  Hawthorne,  et  al.,^'  Hottel,^  Karlovitz,  et  al . ^ Wohl/^ 
and  others.  That  brief  period  provided  some  provacative  experiments  through 
schlieren  photographs^  and  ion  probe  measurements^  that  demonstrated  the 
presence  of  concentrated  reaction  zones  in  turbulent  flames.  Looking  back 

on  some  of  these  experiments,  e.g.,  Hottel,^  Zukoski,^0^  we  find  clear 
evidence  of  the  large  structure  of  inhomogeneous  turbulence  that,  until  the 
recent  work  of  Brown  and  Roshko,^^  escaped  the  attention  of  workers  in 
turbulence  research.  This  early  interest  in  a possible  laminar  flame  struc- 
ture within  turbulent  combustion  zones  stimulated  no  serious  steps  toward  a 

qualitative  formulation  of  these  ideas;  rather,  subsequent  theoretical  devel- 

(12)  (13)  - - (14) 

opments  were  guided  by  the  Reynolds,  Taylor,  Karman,  tradition 

and,  as  a consequence,  treated  the  local  reaction  rates  in  terms  of  mean 

values  of  products  of  fluctuating  state  quantities.  Each  viewpoint  has  its 
merits;  each  has  its  advantage  in  interpreting  the  effects  of  certain 
portions  of  the  turbulence  spectrum,  a fact  that  was  clearly  noted  by 
Damkohler^)  and  Shchelkin.^) 

It  is  the  aim  of  the  present  work,  which  extends  an  analysis  carried 

(15) 

out  by  the  present  authors v ' with  support  of  the  Air  Force  Weapons  Labor- 
atory, to  develop  an  analytical  formulation  to  the  concept  of  a turbulent 
flame  consisting  of  coherent  laminar  flame  elements,  where  by  coherent,  it 

is  implied  that  a local  laminar  flame  element  retains  its  identity  although 
it  is  severely  distorted  and  strained  by  the  turbulent  motions.  More 
specifically,  the  thickness  of  the  diffusion  flame  is  assumed  small  in  com- 
parison with  the  wave  length  of  prominent  disturbances  in  the  gas.  Then 

the  local  flame  element  is  stretched  and  distorted  by  the  local  gas  rotation 
and  strain  rate,  but  the  diffusion  flame  structure  is  affected  only  by  the 
strain  rate  in  its  own  plane. 
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According  to  this  physical  picture,  the  rate  at  which  reactants  are 
consumed,  in  volume  of  dimensions  small  compared  with  the  physical  problem 
but  large  with  respect  to  laminar  flame  thickness,  is  the  product  of  flame 
surface  area  and  reacting  consumption  rate  per  unit  flame  surface.  A 
simplified  illustration  of  this  process  is  shown  in  Figure  1 for  a mixing 

zone.  The  straining  associated  with  the  diffusive  character  of  turbulence 
not  only  leads  to  the  growth  of  flame  surface  area  but  is  essential  in 
determining  the  reactant  consumption  per  unit  of  flame  surface.  A greater 
reaction  rate  is  associated  with  a "stretching"  diffusion  flame  because 
the  straining  in  the  plane  of  the  flame  induces  a flow  of  gas  toward  the 
reaction  line  steepening  the  concentration  gradients  and  thus  augmenting 
the  diffusion  of  reactants  to  the  flame. 


Figure  1.  Schematic  Drawing  of  Flame  Surface  Stretching 
and  Shortening. 

As  the  straining  motion  continues  to  increase  the  flame  surface,  regions 
develop  where  flame  fronts  move  so  close  to  one  another  that  they  interact, 
Figure  1.  In  time,  the  two  neighboring  flame  elements  move  together,  consume 
the  intervening  reactant  and  destroy  these  two  elements  of  flame  surface. 

This  is  the  flame  shortening  mechanism  that  eventually  balances  the  extension 
caused  by  straining,  and  it  is  an  essential  feature  that  was  notably  absent 
from  earlier  intuitive  discussions  of  the  flame  element  model  of  turbulent 


combustion.  It  is  necessary  that  in  the  combustion  of  unmixed  reactants, 
the  competition  between  i)  flame  surface  growth  by  straining  and  ii)  flame 
surface  reduction  by  mutual  annihilation  of  reaction  fronts,  establishes  the 

level  of  reactant  consumption  per  unit  volume. 

A particularly  attractive  feature  of  the  coherent  flame  model  is  the 
effective  separation  of  the  details  of  the  chemistry  from  the  details  of  the 
turbulent  structure.  The  chemistry  appears  only  in  the  plane  laminar  dif- 
fusion flame  subjected  to  strain  rate,  a calculation  in  which  complex  chemistry 
may  be  incorporated  with  relative  ease.  The  turbulent  structure,  on  the  other 

hand,  only  involves  the  heat  release  in  its  response  to  the  corresponding 
density  changes.  Thus,  the  description  of  the  turbulence  and  the  description 
of  chemistry  and  thermodynamics  of  the  reaction  are  coupled  in  a rather 
elementary  way;  when  the  heat  of  reaction  is  low,  the  coupling  is  correspond- 
ingly weak. 


A formulation  of  the  coherent  flame  model  for  the  turbulent  diffusion 
flame  then  requires  i)  a model  for  inhomogeneous  turbulence,  including 
closure  conditions;  ii)  a model  for  flame  surface  distribution  over  the 
turbulent  region,  which  leads  to  a corresponding  reactant  consumption  and 
heat  release;  and  iii)  calculation  or  measurement  of  the  reactant  consump 

tion  rate  for  a laminar  diffusion  flame  that  is  undergoing  strain  in  its 
own  plane. 


In  the  following  sections,  an  analytical  formulation  for  this  model  is 
suggested  and  some  general  properties  are  discussed.  The  equations  are  then 
applied,  for  very  rapid  reaction  rates  and  low  heat  releases,  to  the  problems 
of  the  turbulent  mixing  layer  and  the  turbulent  fuel  jet. 
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2.  THE  COHERENT  FLAME  MODEL 


According  to  the  coherent  flame  model  for  turbulent  diffusion  flames, 
the  flowfield  is  divided  into  two  regions  by  a flame  sheet,  one  of  the 
reactants  existing  on  each  side  of  the  sheet.  The  sheet  becomes  extensively 
distorted  and  dispersed  during  the  turbulent  combustion  process  and  it  is, 
therefore,  appropriate  to  define  a field  variable  z(x.,t)  which  specifies 
the  flame  surface  area  per  unit  volume.  This  quantity  is  denoted  the  flame 
density  and  has  the  dimension  of  a reciprocal  length.  All  of  the  chemical 
reaction  takes  place  within  a region  that  is  small  in  comparison  with  the 
predominant  length  which  describes  the  distortion  of  the  flame  front.  We 

know  that  in  a laminar  diffusion  flame,  such  as  that  represented  by  the 
flame  sheet,  the  combustion  products  remain  within  the  diffusion  layer,  a 
layer  which  under  the  most  elementary  circumstances  grows  as  /t.  The  situ- 
ation is  different  when  the  turbulent  motion  is  continually  extending  the 

flame  surface,  say  at  a linear  strain  rate  t,  and  in  such  a flame,  the 

-1  /2 

diffusion  zones  quickly  stabilize  at  a thickness  proportional  to  e . The 
newly-formed  reaction  products  then  remain  within  this  small  distance  of  the 
flame  sheet  and  are  distributed  over  the  geometric  mixing  zone  by  the  exten- 
sion and  migration  of  the  flame  sheet.  In  the  volume  for  which  the  flame 

density  is  zU^.t),  we  may  similarly  define  mass  fractions  ^(x^t), 

<2(xi  ,t)  of  the  reactants,  the  mass  fraction  k^x^  ,t)  of  the  reaction  pro- 
ducts, and  a mass  fraction  K^(x^,t)  of  an  inert  diluent. 

If  the  process  of  flame  extension  were  to  continue  without  modification, 
the  flame  surface  would  become  dense  in  some  regions  and  the  spacing  between 
surfaces,  which  on  the  average  is  z~\  would  be  of  the  same  order  as  the 
thickness  of  the  flame  diffusion  zones.  In  this  event,  neighboring  flame 

sheets  are  no  longer  isolated;  they  interact  and  quickly  consume  the  inter- 
vening reactant.  ^ This  process  shortens  the  flame  surface  and,  hence, 
reduces  the  flame  surface  density  at  a rate  proportional  to  the  volume  rate 
of  reactant  consumption  by  a unit  area  of  the  flame,  VD»  and  inversely  pro- 
portional to  the  distance  between  the  two  flame  fronts  bounding  the  reactant 
in  question.  Now  the  distance  between  flame  fronts  that  is  occupied  by  a 
particular  species  is  of  the  order  of  the  volume  fraction  of  that  species 
divided  by  the  flame  surface  per  unit  volume.  Then  for  the  species  with 
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volume  fraction  v,-,  the  effective  distance  between  flame  fronts  occupied  by 

i i * 


the  i 


, th 


chemical  species  is  V../E, 


But  is  proves  convenient  to  describe  this 


in  terms  of  mass  fractions  <.  and  the  molecular  weights  M.  of  the  species, 

*i  1 I 

v.  = (k./M.)(Ik./M.)-1,  so  that  the  characteristic  length  associated  with 


.th 


i i 


J J 


the  i chemical  species  is 


K . 


f^-/E  I Kj/Mj 


It  should  be  noted  here  that  when  the  molecular  weights  of  the  reactants 


and  products  are  nearly  equal,  that  is  M = M = M , then  (Xk./M.)  M = 1 

\ u o J J 1 


and  the  above  characteristic  length  becomes  = a form  that  we  shall 

employ  in  later  calculations. 


With  this  result,  we  may  establish  the  physical  order  of  the  rate  of 


.th 


flame  surface  reduction  due  to  the  consumption  of  the  i reactant  as 


VDi 


v . 


per  unit  of  flame  surface,  or 


VDi  e2 


(1) 


l 


per  unit  of  volume  that  contains  a flame  area  e.  We  shall  write  the  flame 
surface  reduction  caused  by  consumption  of  the  ftfel  and  oxidizer  components, 
i * 1 and  i = 2,  as  the  sum  of  terms  of  the  type  given  by  expression  1. 


Now  it  is  a matter  of  formal  calculation  to  obtain  an  equation  for  the 
change  in  area  of  an  element  associated  with  fluid  particles,  as  it  is 
stretched  and  deformed  by  a turbulent  medium.  In  fact,  Batchelor^  ^ has 
done  exactly  this  calculation  in  examining  the  behavior  of  line  segments 
and  surface  elements  in  homogeneous  isotropic  turbulence.  For  our  purposes, 
the  result  states  that  the  Eulerian  time  derivative  of  E is  given  by  the  sum 
of  the  turbulent  diffusion  of  this  quantity  and  the  increase  of  flame  density 
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by  the  local  strain  rate  of  the  mean  motion.  It  is  presumed  that,  for 
inhomogeneous  turbulence,  the  strain  rate  of  the  mean  motion  is  sufficiently 

dominant  that  the  additional  straining  associated  with  turbulent  fluctuations 
may  be  neglected.  For  two  reasons,  this  result  does  not  describe  accurately 

the  behavior  of  a flame  element.  First,  the  flame  front  generally  moves  with 
respect  to  the  fluid  and  thus  does  not  always  contain  the  s^m-a  fluid  particles 
Second,  the  flame  shortening  mechanism,  in  which  adjacent  flame  fronts  consume 
the  intervening  reactant  and  annihilate  each  other,  must  be  accounted  for. 

They  are  two  aspects  of  the  same  differences  which  contrast  the  behavior  of 
the  flame  sheet  with  that  of  a surface  element  consisting  of  the  same  fluid 
elements.  The  first  of  these  will  be  neglected  in  our  formulation  while  the 
second,  and  dominant,  effect  will  be  accounted  for  by  shortening  mechanisms 

discussed  above. 

Then  in  a turbulent  fluid  field  with  mean  velocity  components  IL , scalar 
turbulent  diffusivity  V and  local  scalar  mean  strain  rate  e,  we  postulate 

that  the  flame  surface  density  satisfies  the  partial  differential  equation 


where  the  specific  subscripts  j = 1 and  j = 2 will  denote  the  "fuel"  and 
"oxidizer"  reactants,  respectively. 

A notable  advantage  of  the  coherent  flame  model,  upon  which  we  shall 
elaborate  later,  is  that  the  chemical  reaction  rate  may  be  expressed  in 
terms  of  the  local  flame  surface  density.  Thus,  for  example,  the  mass 
of  species  1 that  is  being  consumed  per  unit  volume  is  just  p Z since 
Voi  is  volume  consumption  rate  of  fuel  by  a unit  flame  area.  Then  by  well 
known  techniques,  we  write  the  continuity  equation  for  the  fuel  constituent 
in  a turbulent  flow 


6 


and  similarly  for  the  oxidizer  component. 


3<2 
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where  a*  is  an  appropriate  turbulent  Schmidt  number.  Because  £ kj  = 1>  the 
remaining  species  conservation  equation  is  unnecessary  so  long  as  the  inert 
diluent  may  be  grouped  with  other  species  or  is  absent.  It  is  assumed,  of 
course,  that  the  reactant  consumption  rates  V^,  V 2 for  a unit  flame  area 

are  available  in  terms  of  the  local  variables  and  transport  properties. 

The  reactant  consumption  rate  may,  in  the  case  of  the  complex  chemistry, 
require  a numerical  calculation  or,  in  the  case  of  rapid  kinetics  which  we 
shall  examine  subsequently  in  more  detail,  it  may  be  expressed  analytically 
in  a simple  closed  form.  Under  any  circumstances,  the  quantity  of  importance, 
Vp-j , is  obtained  from  a one-dimensional  calculation  that  is  independent  of 

the  detailed  solution  of  the  partial  differential  equations  describing  the 
turbulent  flame. 

To  complete  the  formulation,  a model  for  inhomogeneous  turbulence  is 
required  which  yields,  in  addition  to  the  velocity  field,  the  scalar  dif- 
fusivity  and  the  scalar  straining  rate  e.  The  model  that  is  the  most  closely 
related  is  that  suggested  by  Saffman^)  wbich  utilizes  the  specific  vorticity 
H and  the  specific  turbulent  kinetic  energy  as  independent  variables, 
employs  a scalar  diffusivity  and  recognizes  the  role  of  the  local  magnitude 
of  the  straining  motion  in  the  mechanism  for  producing  both  vorticity  and 
energy.  Then  for  an  incompressible  medium,  implying  that  the  heat  of  reac- 
tion is  small,  this  model  gives 
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The  first  two  equations  are  of  the  standard  form,  and  it  is  assumed  that 
the  Reynolds  stresses  — as  well  as  other  diffusion  processes  --  may  be 
described  by  a scalar  turbulent  diffusivity  V.  It  is  a familiar  dimensional 
argument  that  the  turbulent  diffusivity,  which  has  the  dimension  of  a veloc- 
ity times  a length,  may  be  constructed  from  local  field  quantities;  in  fact, 

this  is  the  essence  of  Karman’s  original  concept  of  mechanical  similitude. 

2 2 

In  the  model  we  use,  Saffman  chooses  the  turbulent  energy  density  e («.  ft  ) 
and  the  magnitude  oj  of  the  vorticity  (1/t)  as  the  local  physical  quantities. 
Both  turbulent  kinetic  energy  and  vorticity  have  physically  understood 
transport  laws,  and  they  are  related  through  the  association  of  the  turbulent 

velocities  with  the  vorticity  that  induces  them.  Then  take 


V = e/w  (9; 

which  not  only  defines  the  turbulent  diffusivity,  but  establishes  a magni- 
tude relationship  between  e and  w. 

The  transport  Equation  7 for  the  turbulent  energy  has  a physically 
recognizable  form;  the  three  terms  on  the  right-hand  side  are,  in  order; 
diffusion,  production,  and  dissipation.  The  turbulent  diffusion  term  is 
clear  and  the  production  term  is  patterned  after  the  one  which  occurs 

naturally  In  turbulence  theory  which  is  proportional  to  the  Reynolds  stresses 

acting  on  the  mean  strain  rates.  Here,  e is  taken  as  a measure  of  the 

Reynolds  stress,  which  j (9U^/5Xj  + aUj / 3x ^ ) is  the  strain  rate  of  the  mean 

flow,  and  we  have  denoted  it  S...  The  last  term,  the  dissipation,  must  be 

ij 

constructed  from  the  product  of  three  fluctuating  velocities,  and  Saffman 
has  chosen  this  to  be  the  product  of  vorticity  and  turbulent  energy.  The 
vorticity  equation  8 is  reminescent  of  Helmholtz  vortex  equations  and  may 
be  intuitati vely  argued  from  them.  Again,  the  turbulent  diffusion  of 


vorticity  Is  described  by  the  diffuslvity  V.  The  second  term  represents  the 
production  of  vorticity  by  vortex  stretching,  the  stretching  being  given  In 
Sflff man's  model  by  the  magnitude  of  the  velocity  gradient  tensor.  It  should 
be  noted  that,  the  mean  velocity  derivatives  enter  differently  in  the  energy 
and  vorticity  production  terms,  Saffman,^®)  and  it  is  only  for  certain 

flowfields  that  the  two  become  indistinguishable. 

This  set  of  four  equations,  5 through  8,  results  from  a complete  model 
in  the  sense  that  the  constants  involved  are  universal,  closure  being 
achieved  by  the  choice  of  turbulent  diffuslvity  V s e/w.  In  the  original 
formulation  of  the  coherent  flame  model  presented  in  Reference  15,  closure 
was  achieved  by  identifying  the  turbulent  diffusivity  p as  the  product  of 

a characteristic  velocity  and  a characteristic  length  defined  locally  or 
globally.  The  fact  that  this  choice  must  be  made  with  some  physical  insight, 
and  differently  for  each  type  of  flow,  shows  that  the  problem  is  not  des- 
cribed completely  by  the  differential  equations;  such  a model  is  sometimes 

denoted  incomplete,  implying  that,  when  coupled  with  the  coherent  flame 
model,  the  coefficients  a,  a*,  etc.  are  not  universal  but  may  require  adjust- 
ment from  one  flow  type  to  another. 
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3.  REACTANT  CONSUMPTION  RATE 

The  reactant  consumption  rates  enter  into  the  problem  in  the  form  VD-| 
and  Vq2>  the  volume  of  fuel  and  oxidizer  consumed  ptr  unit  flame  area, 
respectively,  and  play  roles  both  in  the  species  conservation  relations  and 
in  the  equations  describing  the  flame  surface  density.  These  quantities 
are  assumed  to  be  determined  locally  by  the  flame  structure  and  to  depend 
only  upon  local  quantities;  in  particular,  for  the  diffusion  flame,  they 
are  determined  by  fuel  and  oxidizer  concentrations  arid  a local  fluid  mech- 
anical property.  As  shall  be  illustrated,  they  may  be  determined  analytic- 
ally, or  by  numerical  calculations  when  the  kinetics  are  complicated  and 
have  an  essential  effect.  The  important  point  to  keep  in  mind  is  that  the 
entire  flame  structure  and  chemical  kinetics  are  coupled  with  the  field 
analysis  rather  weakly,  so  that  the  consideration  of  complex  kinetics 
complicates  only  the  local  flame  structure  and  not  the  formulation  of  the 

flowfield  and  the  flame  distribution. 

As  a first  illustration,  consider  the  diffusion  flame  with  rapid 
kinetics;  in  this  approximation  the  reaction  takes  place  at  a surface  of 
infinitesimal  thickness.  Utilizing  the  coordinates  x and  y to  signify 
distances  parallel  to  and  normal  to  the  flame  surface,  supposed  to  lie 
along  y = 0,  the  fuel  and  oxidizer  concentrations  satisfy  the  equations 


where 
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when  the  heat  of  reaction  is  negligible  and  the  gas  density  p is  a constant. 
For  the  classical  diffusion  flame,  the  field  is  independent  of  x,  but  the 


solution  is  time  dependent  Thoti  c'in‘"Q  8u 
appear,  is  not  generally  zero.  Then  introducing  the  variable 


Then,  since  ly.  = 0,  v = v(t)  which,  as  will 

3X 
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03) 


and  taking 


v(t)  = W 


04) 


the  species  conservation  equations  are  reduced  to  a similarity  form  and 
become 


d2tcj 

S7 


(H 


dx.j 

df 
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(15) 


where  i = 1,  2 for  fuel  or  oxidizer,  respectively.  This  pair  of  differ- 
ential equations  is  required  to  satisfy  the  conditions  that  the  fuel  and 
oxidizer  mass  fractions  take  on  the  values  k^®)  and  k2(-°°)  at  y * +» 
and  y = -°°,  respectively  and  that  the  mass  flux  to  the  diffusion  flame, 

y = 0,  supplies  fuel  and  oxidizer  in  the  stoichiometric  ratio.  This  latter 
condition  is  explicitly 


PD  fly  (0»t) 

- pD  ^ (O.t) 


06) 


where  f is  the  known,  constant  stoichiometric  fuel-oxidizer  ratio.  It  is 
not  difficult  to  show  that 


and 


erf  (|  - w)  + erf  (W) 
K1  = Kl("}  rVerf  7w) 


erf  (-  | + w)  - erf  (W) 
i . erf  (W) 


(17) 


(18) 


satisfying  the  differential  Equation  15  and  the  boundary  conditions  at 
y - ±«.  The  stoichiometry  condition  (Equation  16)  then  determines  the 
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characteristic  value  W,  and  it  is  a matter  of  calculation  to  show  that  this 
gives  the  result 

<1 (°°)  [l  - erf  (W) 

^(-“O  ^1  + erf  (W) 

Now  k-|  (~)/k2 (-00)  is  the  imposed  fuel-oxidizer  ratio  of  the  problem,  aid  the 
quotient  of  this  with  the  stoichiometric  fuel-oxidizer  ratio 


K.j  (00)/k2(-co) 

f 


E <f> 


(20) 


is  frequently  called  the  equivalence  ratio.  Thus,  the  value  of  W becomes 


This  quantity  defines  the  value  of  the  transverse  gas  velocity 


(21) 


v(t)  - ys  erf-1  (f^4)  (22) 


which  is  required  to  keep  the  flame  stationary  at  the  x-axis.  With  this 
solution,  the  values  of  the  reactant  volume  consumption  rates  may  be  cal- 


culated as 


and 


VD1=k1(")  ^(*2^-) 


-w2 


VD2  = k2(’°o) 


-W2 


(23) 


(24) 


This  result  exhibits  the  intuitively  clear  result  that  the  consumption 
rates  decrease  with  increasing  time,  because  the  diffusion  layers  that 
supply  the  reactants  grow  thicker  with  time.  The  equivalence  ratio,  which 
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determines  the  value  of  W,  is  known  because  (°°)  and  ^(-w)  are  equal  to 
the  values  remote  from  the  turbulent  flame  since  the  diffusion  zone  thick- 
ness is  assumed  small  in  comparison  with  flame  front  spacing.  At  any  point 

within  the  turbulent  flame,  therefore,  the  reactant  consumption  rates  are 
known  in  terms  of  the  time  t elapsed  since  the  formation  of  the  flame. 

In  disucssing  the  turbulent  flame  structure  earlier,  it  was  emphasized 
that  the  turbulent  motions  tend  to  extend  the  flame  surface  and  that  the 
significant  part  of  this  extension  consists  of  strain  rate  in  the  plane  of 
the  flame.  Now,  if  the  flame  is  aligned  with  the  x-axis  and  the  remaining 
axes  chosen  so  that  straining  rate  is  along  the  x-axis,  the  resulting 
strained  diffusion  flame  is  also  directly  soluble.  In  the  particular 
instance  when  the  straining  rate  in  the  fluid  is  large,  i.e.,  where 
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- M 

“ 3X 


(25) 


is  large  and  constant.  Equations  10  through  14  take  the  form 


and 


(26) 


(27) 


If  we  solve  the  problem  as  previously  posed,  the  solution  has  a brief 
transient  followed  by  a steady  solution  which  we  may  obtain  be  neglecting 
3K.j/9t  and  taking  v and  to  be  functions  of  y only.  Then  introducing 
the  variable 


? = 

/D/2e 


v = -ey  + /2cD  W 

with  W unknown,  the  differential  equation  26  becomes 
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formally  identical  with  Equation  15.  The  conditions  for  y ± »,  as  well  as 
stoichiometry  condition  at  the  flame  front,  are  also  identical  with  those 
of  the  time-dependent  flame,  and,  hence,  the  solution  for  the  strained 
flame  is  likewise  given  by  Equations  17  and  18  together  with  W evaluated 
by  Equation  21,  but  with  £ as  defined  in  Equation  28.  The  corresponding 

volume  consumption  rates  of  reactants  are  easily  calculated 


(3D 


Thus,  so  far  as  the  local  flame  structure  is  concerned,  the  features 
relevant  to  the  turbulent  flame  calculation  are  determined  by  the  local 
strain  rate.  Because  of  the  simplicity  of  these  analytical  results,  they 
will  be  employed  in  detailed  calculations  for  the  mixing  zone  and  the  fuel 
jet. 

In  many  physical  problems,  the  detailed  distribution  of  chemical 
species  plays  a significant  role,  and,  consequently,  the  simplification  of 
the  preceding  calculation  cannot  be  employed.  But  because  of  the  effective 
separation  of  chemical  aspects  and  turbulence,  numerical  calculations  of 
even  very  complex  flame  structures  is  possible.  To  illustrate  such  a com- 
putation for  a complex  chemical  reaction,  we  consider  the  reaction  between 
hydrogen  and  atomic  fluorine  yielding  hydrogen  fluoride  in  the  vibrational 
states  v:  . 

H2  + F + HF(v)  + H v = 0,  1 , 2,  and  3. 


Because  the  vibrational  population  distribution  produced  by  the  reaction 
is  partially  inverted,  the  hydrogen  fluoride  can  be  caused  to  lase  (in  the 


'.m 
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infrared)  in  a properly-designed  laser  cavity.  In  fact,  H2  - F and  D2  - F 
chemical  lasers  are  among  the  most  efficient.  For  this  reason,  the  important 
rate  coefficients  have  been  determined  to  considerable  accuracy  and  for  the 
same  reason,  there  is  great  interest  in  the  turbulent  combustion  of  these 

reactants.  It  is  natural,  therefore,  to  choose  the  H2  - F flame  for  the 
illustration. 


It  is  the  chemical  laser  application  as  well  that  necessitates  numerical 
calculation  of  the  flame  structure.  To  determine  laser  performance,  it  is 
not  sufficient  to  know  the  and  F consumption  rates.  In  addition,  the 
production  of  HF  in  the  various  vibrational  levels  is  needed  as  is  the  sub- 
sequent depopulation  rate  of  these  levels  by  inter-molecular  collisions. 

The  accurate  treatment  of  these  processes  and  of  the  inter-diffusion  of  the 
many  species  involved  is  possible  only  in  a numerical  calculation. 


For  reasons  irrelevant  to  the  present  discussion,  the  computations 
were  carried  out  utilizing  parallel  streams  of  reactants  that  cause  the 
development  of  the  flame  along  the  axis  of  the  mixing  zone.  When  the  two 
streams  have  the  same  velocity  U,  the  successive  values  of  x can  be  con- 
sidered as  successive  intervals  in  time  defined  by  t = x/U.  The  results 
then  become  comparable  to  the  time  dependent  diffusion  flame  with  rapid 
kinetics  discussed  above.  Correspondingly,  when  a pressure  gradient  is 
imposed  upon  the  parallel  stream  to  accelerate  them  in  their  direction  of 
motion,  the  flame  eventually  reaches  a structure  that  is  independent  of 
x and  directly  comparable  with  the  steady,  strained  solution  for  rapid 
kinetics,  also  discussed  above. 


Consider  then  parallel,  accelerating  streams  of  nitrogen,  one  carrying 
the  hydrogen  and  the  other  atomic  fluorine.  See  Figure  2.  At  the  initial 
station,  x = 0,  the  pressure  in  both  streams  is  5 torr,  the  temperature  is 
500°K,  and  the  velocity  10  cm/sec.  In  the  upper  stream,  the  fluorine  and 
nitrogen  mass  fractions  are  0.077  and  0.923,  respectively;  in  the  lower, 

the  hydrogen  and  nitrogen  mass  fractions  are  0.011  and  0.989,  respectively. 


A constant  strain  is  imposed  on  the  flow  by  causing  both  streams  to 
• 4 

accelerate  uniformly  to  10  cm/sec  in  a distance  of  6 inches.  In  the 


process,  the  pressure  and  temperature  drop  to  4.67  torr  and  490°K,  respec- 
tively. 
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Figure  2.  Strained  Laminar  Mixing  Layer 


The  chemical  reactions  that  are  considered  are  described  in  detail  in 
Reference  20.  They  may  be  summarized  as  follows.  The  reaction  + F -v 

HF ( v ) + H is  unidirectional  and  produces  vibrational  populations  in  the 
following  fractions:  v = 0,  0.056;  v = 1,  0.111;  v = 2,  0.555;  and 

v = 3,  0.278.  The  excited  states  are  deactivated  by  collision  with  the 
other  species  present  at  prescribed  efficiencies.  In  addition,  there  are 

V-V  exchange  collisions  among  the  various  levels  of  HF  and  with  molecular 
hydrogen  and  nitrogen.  The  rate  coefficients  for  all  these  molecular  pro- 
cesses are  given  in  Reference  20. 

The  multi -component  diffusion  coefficients  are  computed  from  temper- 
ature dependent  binary  coefficients  under  the  assumption  that  the  Lewis- 
Semenov  number  was  constant.  The  temperature  dependence  of  the  species 
viscosities  is  accounted  for  by  empirical  fits  to  experimental  data. 

The  system  of  equations  describing  this  diffusing  and  chemically- 
reacting system  was  solved  by  a computer  program  that  is  a modification  of 
the  Blottner  boundary  layer  code  described  in  Reference  21.  The  original 
program  has  been  modified  so  that  it  solves  the  laminar  mixing  layer 
equations  and  has  subsequently  been  used  in  several  investigations  of  the 
of  the  constant  pressure  hydrogen-fluorine  flame. A further  modifica- 
tion was  made  to  allow  the  treatment  of  the  strained  flame  in  the  present 
study. 


*The  authors  thank  F.  E.  Fendell  and  D.  Haflinger  for  making  this  modification 
and  for  obtaining  the  computer  results  that  are  presented. 
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Two  sets  of  calculations  were  made:  one  with  the  complete  chemical 

system  outlined  above  and  a second  in  which  the  HF  was  taken  to  be  produced 

in  a single  (ground)  state.  This  latter  simplification  was  made  so  that  an 
investigation  of  the  effect  of  rate  coefficient  changes  could  be  made  more 
economically.  These  results  are  discussed  first. 

Figures  3 through  5 show  the  Hg,  F and  HF(0)  y distributions  at  several 
downstream  locations;  that  for  Figure  3 being  close  to  the  initial  stations. 
We  observe  first  the  evolution  of  the  profiles  into  the  fixed  shape  that 
the  simple  analysis,  outlined  at  the  beginning  of  this  section,  predicts 
for  the  strained  flame.  The  profiles  at  x = 6.0  inches  are  virtually  the 
same  as  those  at  x = 0.126  inch.  This  asymtotic  state  is  attained  even 
though  neither  the  density  nor  the  velocity  are  constant  in  the  reaction 
zone.  The  density  is  reduced  as  the  heat  of  reaction  increases  the  temper- 
ature and  then  this  lighter  gas  accelerates  more  than  the  free  stream  under 
the  imposed  pressure  gradients. 

The  next  point  of  interest  in  these  results  has  to  do  with  the  chemical 
reaction  rate  relative  to  diffusion  rate.  We  note  that  even  for  x locations 
as  close  as  0.126  inch  to  the  origin,  there  is  little  interpenetration  of 

hydrogen  and  fluorine.  This  observation  suggests  that  the  assumption  made 
in  the  earlier  analysis  that  the  reaction  rate  is  infinite  may  be  a good 
approximation  for  the  present  condition.  Additional  evidence  that  this  is 
the  case  was  obtained  by  repeating  the  computation  with  the  reaction  rate 

coefficient  increased  by  a factor  ten.  Comparison  of  Figure  6,  given  by 
the  computation  with  Figure  5,  shows  only  a minor  increase  in  the  amount 
of  HF  that  has  been  formed,  a result  implying  that  diffusion  to  the  flame 
surface  is  the  rate-limiting  process. 

Vibrational  state  distributions  for  the  run  with  the  full  hydrogen- 

fluorine  chemistry  are  given  in  Figure  7 through  10.  We  see  that  the 

vibrational  population  total  inversion  which  exists  near  the  origin, 

x = 10"®  inches,  has,  under  the  action  of  the  deactivating  collisions, 

« -5 

nearly  disappeared  at  x = 4.4  • 10  inches.  (The  inflection  in  the  HF(v) 
distribution  at  the  very  low  concentrations  comes  from  computation  errors 
and  may  be  ignored.) 
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Distributions  in  Strained  Figure  6.  H2,  F,  and  HF  Distributions  in  Strained 

lame  at  x = 6.0  Inches.  Laminar  H2-F  Flame  at  x = 6.0  Inches  with 

Reaction  Rate  Coefficients  Increased  by 
Factor  10. 


Ground  and  Excited  HF  Vibrational  State  Figure  10.  Ground  and  Excited  HF  Vibrational  State 
Distributions  in  Strained  Laminar  H2-F  Distributions  in  Strained  Laminar  H2-F 
Flame  at  x = 7.29  * 10-2  Inches.  Flame  at  x = 0.28  Inches. 


Comparison  of  Figures  9 and  10  indicates  that,  in  this  case,  also  the 
imposed  constant  strain  rate  has  formed  a flame  whose  width  is  constant  and 

within  which  the  state  distributions  are  nearly  constant.  This  result  is 
perhaps  somewhat  unexpected  in  that  the  deactivation  processes  are  distri- 
buted throughout  the  diffusion  zone  and,  thus,  the  conditions  of  the  chemi- 
cal reaction  are  quite  different  from  those  in  the  simple  analysis  that  lead 

to  the  constant  state. 

These  results  would  be  incorporated  in  the  coherent  flame  model  by 
determining  the  fuel  and  oxidizer  consumption  rates  from  the  numerical 
solution.  The  state  of  the  fluid  within  the  turbulent  flame  zone  would  then 
be  fixed  by  knowledge  of  the  flame  surface  per  unit  volume,  determined  from 
the  model,  and  the  distribution  within  the  flame.  For  all  axial  locations 

except  those  very  close  to  the  origin,  the  constant  state  strained  distri- 
butions would  be  the  appropriate  ones.  It  appears,  in  fact,  that  in  most 
applications,  the  strain-dominated  flame  solutions  would  be  used. 

In  both  the  analysis  and  the  numerical  calculations  just  described, 
the  assumption  is  made  that  no  products  are  present  when  the  flame  begins. 

It  is  clear,  however,  that  in  those  regions  of  the  turbulent  zone  where 
flame  shortening  occurs,  products  are  left  behind  which  may  influence  the 
subsequent  reaction.  They  may,  for  instance,  separate  volumes  of  fresh 
reactant  and  thereby  reduce  the  reaction  rate.  The  model  has  not  yet  been 

modified  to  account  for  such  effects,  but  some  indication  of  their  signi- 
ficance may  be  obtained  from  the  following  idealization.  Consider,  as  an 
initial  condition,  the  situation  sketched  in  Figure  11,  meant  to  represent 


concentration 


Figure  11.  Initial  Distribution  of  Reactants  Separated  by  Product 
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fuel  ami  oxidizer  separated  by  an  Inert  product  layer  2 xQ  thick.  With 
the  simplifying  assumptions  that  the  stoichiometric  ratio  is  unity,  the 
fuel  and  oxidizer  concentrations  are  equal,  the  rate  is  infinite,  and  the 

density  Is  constant,  the  reaction  is  confined  to  the  plane  x » 0 and  the 
reactant  concentrations  there  are  held  at  zero  for  t > 0.  The  fuel  con- 
centration distribution  satisfying  the  time-dependent  diffusion  equation 
with  these  boundary  and  Initial  conditions  is: 
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which  can  be  written 

x +x/2fUf 


x -x/2v/Uf 
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A similar  expression  applies,  of  course,  to  the  oxidizer.  From  these 
solutions,  the  reactant  flux  to  the  flame  sheet  is  found  to  be 
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from  which  it  is  seen  that  the  consumption  rates  differ  from  those  of  the 
classical  case  by  the  factor  exp  (-x  /4Dt).  This  result  makes  it  clear 
that  when  large  masses  of  product  are  interposed  between  the  reactants, 
the  diffusion  flame  structure  may  be  strongly  altered. 
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4.  REACTION  IN  A TURBULENT  MIXING  ZONE 


Consider  the  turbulent  mixing  and  reaction  of  two  parallel  streams, 
the  upper  with  undisturbed  velocity  Up  fuel  mass  fraction  k^(«)  and  diluent 
mass  fraction  k^(«)  = 1 - k-j(“),  the  lower  with  undisturbed  velocity  U2» 
oxidizer  mass  fraction  ^(-w)  and  diluent  mass  fraction  = 1 - 

The  problem  is  two-dimensional,  steady,  and  will  be  treated  under  the 
boundary  layer  approximation.  If  we  assume  further  that  the  molecular 
weights  are  nearly  equal  and  that  the  molecularity  of  the  reaction  and 
heat  evolved  are  small  enough  to  produce  negligible  volumetric  changes, 
then  the  coherent  flame  leads  to  the  following  formulation  of  the  problem 
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In  this  problem,  the  magnitude  of  the  mean  rate-of-strain  tensor  and  the 
magnitude  of  the  velocity  gradient  tensor,  as  employed  by  Saffman,  have 
the  same  value  and  appe?  in  Equations  35  and  36.  The  turbulent  exchange 
ratios,  which  relate  the  turbulent  diffusion  of  energy  or  vorticity  to  the 
turbulent  diffusion  of  momentum,  have  each  been  taken  as  1/2,  values  found 
acceptable  in  Saffman's  calculations. 

The  chemical  reaction  portion  of  the  model  will  be  formulated  here 
utilizing  the  flame  structure  based  upon  very  fast  kinetics  and  dependence 
upon  the  local  strain  rate.  According  to  this  analysis,  given  in  Section  3, 
the  local  volumetric  consumption  rates  of  fuel  and  oxidizer  per  unit  flame 
area  are  i^H  ■—  g(+)  aU/ay | and  k2(-«)  g(<f>)  /D|  au/ay ( where  again  ( aU/ay | 
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appears  as  the  magnitude  of  the  strain-rate  tensor,  a quantity  denoted  e in 
Section  3 and  g(<j>)  = (1  + <}>)  exp  (-W2).  Then  the  conservation  equations 
for  the  fuel  and  oxidizer  are 


'i<-)}g<*>  AifT  = 

k2(-«>)  g(4> ) D | |y I t. 


(37) 


(38) 


For  the  mass  fractions  of  the  product  and  inert  diluent,  we  have  similarly 


J!3  . ,,  _ 3__  e ^3 

3x  v 3y  ' 3y  a)  3y 


K-|((j))  + 


(39) 


U 


9X 


+ V 


(40) 


Because  of  our  assumptions  regarding  the  nearly  identical  molecular  weights 
of  the  gasses,  the  inert  diluent  is  indistinguishable  from  other  components 
and,  hence,  may  be  grouped  with  any  of  the  other  constituents.  The  only 
alteration  will  be  to  define 
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and  Kg  so  as  to  satisfy  the  relation  ^ ^ = 1.  This  general  rela- 

tion permits  us  to  ignore  Equation  39  and  to  rewrite  37  and  38 
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Finally,  the  flame  surface  density  equation  is 


U 3£  4.  y 3E 
3X  3y 


(*  &)+  a 
\<*>  3 y/ 


- X igU)  l/oi^T  ^ 
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Here  we  have  used  the  same  universal  constant  a in  representing  the  flame 
surface  stretching  as  was  used  in  representing  the  vortex  stretching  in 
Equation  36.  Note  further  that  we  have  introduced  only  one  additional 
universal  constant,  A,  in  the  flame  shortening  terms.  The  quadratic  nature 
of  the  shortening  terms  allows  any  universal  multiplier  to  be  abosrbed  in 
the  definition  of  £ itself;  hence,  only  one  additional  constant  is  necessary. 

Now  because  of  our  assumption  of  small  density  changes,  the  fluid 
mechanical  problem  is  uncoupled  from  the  flame  structure,  and  a similarity 
solution  for  the  turbulent  mixing  zone  may  be  obtained  in  the  same  manner 
as  given  by  Saffman.  Then  introducing  £ = x,  n = y/x  and  defining  a stream 
function 

* = U*  5 F(n)  (44) 


where 


U*  = (U-]  + Ug) 


the  velocity  components  become 


U * U*  F'  U)  {46, 

V - U*  (n  F*  - F) 

Through  dimensional  considerations,  it  is  appropriate  to  write  the  vorticity 
and  the  energy  densities  as 
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(d  = U*  j a (n) 


(47) 


e 


(n) 


(48) 


where  n and  E are  dimensionless  functions  of  the  similarity  variable  alone. 
With  this  representation,  the  turbulent  diffusivity  is 


(49) 


a result  which  allows  the  momentum  Equation  34  to  be  written  in  similarity 
form 


F F"  + <L  (L  F"\  = 0 

dn  ) 


(50) 


With  the  similarity  formalism  already  introduced,  the  vorticity  and 
energy  density  equations.  Equations  35  and  36,  reduce  to  the  following 
ordinary  differential  equations 


and 


(Ffl)  = ^ (e^‘  ) + a F" 


6 ft3 


-F  E' 


+ a*  F"E 


3*  a E 


(51) 


(52) 


Equations  50,  51,  and  52  correspond  to  those  given  by  Saffman  for  the 
mixing  zone  problem. 

Turing  now  to  the  equations  related  to  the  flame  and  combustion 
problem,  it  is  not  obvious  that  they  permit  a corresponding  similarity 
representation.  The  peculiar  term  that  occurs  in  Equations  41,  42,  and  43 
is  /D|aU/3y|  which  we  may  write  in  terms  of  the  similarity  variables  as 
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Volf ! - V°i-If  <U*  F')|  - U,  V4  |F"(n)|  . (53) 

Then  the  formal  representation  of  the  flame  structure.  Equations  41  through 
43,  in  the  similarity  variables  of  the  fluid  mechanical  problem  gives 
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(55) 


(56) 


By  inspection  of  these  results,  it  appears  that  for  similarity  ie  and 
are  independent  of  5 and  that  £ behaves  as  Thus,  we  have 


K-|(n):  <i  (00)  = 1 , K^(-«>)  = 0 


^(h)-'  Kg (°°)  “ 0 > k2(-®)  - 1 


(57) 

(58) 


and  the  flame  density  is  expressed  as 


£ U.n)  = 


L(n) 


9U)  flm: 


(59) 


Note  the  appearance  in  this  expression  and  in  Equations  54  through  55 
of  the  dimensionless  streamwise  variable 


28 


C = 


(60) 


which  bears  the  same  relation  to  the  Reynolds  number  as  the  molecular  dif- 
fusivity  D does  to  the  kinematic  viscosity.  Note,  then,  that  the  quantity 
/DS/U*  has  the  dimension  of  a length  and  gives  the  flame  surface  density  its 
dimension  of  a reciprocal  length,  the  variable  L ( n ) being  the  dimensionless. 
Then,  we  use  the  new  density  variable 


L(n) : L(“)  = L(-“)  = 0 (61 ) 

With  these  definitions.  Equations  54  through  56  reduce  to  ordinary 
differential  equations  in  k (n),  ^(n),  and  L(n): 


In  contrast  to  some  formulations  of  turbulent  mixing  problems,  the 
present  one,  because  of  its  non-linear  turbulent  diffusivity,  may  produce 

a sharp  boundary  between  the  turbulent  zone  and  the  far  field.  This  not 
unfamiliar  situation  was  discussed  in  some  detail  by  Saffman^8^  in  connection 
with  the  fluid  dynamic  solutions  of  several  problems.  The  relevant  point 
here  is  that,  with  the  non-linear  diffusivity  employed,  the  edges  in  this 

case  are  sharp,  lying  at  constant  values  of  n,  these  values  emerging  as 
characteristic  of  the  problem.  With  the  value  of  1/2  chosen  as  the 
exchange  ratio  in  Equations  35  and  36,  the  appropriate  boundary  values 
for  F,  E,  and  n are 
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U7 

F(ni ) = u; 

u2 

F(n2)  - ^ 

fi(n-j ) = fi(n2)  = 0 

E^)  = E‘  (n j ) = E(n2)  = E'(n2)  = 0 


(65) 


Subjecting  the  flame  variables  , k2>  and  L to  a corresponding  analysis, 
we  find  for  the  equations  — as  formulated  — 


K1 (n1 ^ " K1 = 0 

^2(nl)  = ^2 ( n2  ^ = 0 
L(n-|)  = L(n2)  = 0 


(66) 


Clearly,  there  exists  a trivial  solution  L(n)  = 0,  and  in  this  case,  tc 


1 


and  k2  simply  become  constituents  that  mix  without  chemical  reaction.  In 
fact,  K-|(n)  and  ic2(n)  are  linear  functions  of  F'(n),  but  this  is  not  the 
solution  we  seek. 


Some  properties  of  the  solution  can  be  obtained  quite  easily,  for 
example  the  total  fuel  consumed  per  unit  length  of  the  mixing  zone.  This 
follows  from  integration  of  either  Equation  41  or  Equation  62  across  the 
mixing  layer.  Working  with  Equation  41,  we  find 


yi  yi  • yl  

j'  ax  (U  Ki)  dy  + J'  Jy  (v  Ki)  dy  = ” J'  ^9(4)  V°i  ayl  E’  dy 

H y2  y2 


which,  with  appropriate  substitution  for  the  variables,  leads  to 
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dx 
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I 

dy  - U)n1  - 1 U*  y 4TT  L dn 


(67) 


This  expresses  the  change  with  respect  to  x of  fuel  flux  in  the  mixing 

zone,  given  by  the  integral  on  the  left,  as  the  difference  between  the 

flux  of  fuel  in  at  the  upper  boundary  and  the  fuel  consumed  by  the  chemi- 
cal reaction.  The  striking  result  is  that  this  rate  is  independent  of 

distance  along  the  mixing  zone,  in  spite  of  the  fact  that  all  other 

physical  quantities  have  some  x-dependence. 


Calculations  by  Milinazzo  and  Saffman^8)  have  provided  a good  basis 
for  the  formulation  of  an  integral  solution  to  the  turbulent  mixing  zone 
which  we  carry  through  the  flame  zone  problem  as  well.  For  this  purpose, 
we  approximate  the  stream  functions 
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U*  n2 


n + 


U1  U2 

, tc-n;  2 
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(68) 


fl(n)  = (n-|  " n)  (n  - n 2} 


(69) 


E(n)  = E0  (n-|  - n)2  (n  - r^)2 


(70) 


in  which  n-j , n2>  ^0»  and  Eq  appear  as  constants  (functions  of  the  problem 
parameters)  to  be  determined  by  the  integral  relations.  Using  integrals 
of  the  energy  and  vorticity  equations,  together  with  integrals  from  0-n 
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and  from  03-0  of  the  momentum  equation,  the  problem  reduces  to  an  algebraic 
solution.  The  technique  is  familiar  and  straightforward,  although  alge- 
braically tedious  and  the  details  will  be  omitted.  For  small  values  of 
we  find  that 


(U]  - u2)/u, 
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EJ  , — 1 — } = Q0  = 153.64 


•°U  (u,  +u2); 


(71) 


The  values  of  n^  for  specific  values  of  U2/U-]  correspond  very  well  with 
numerical  results  given  by  Milinazzo  and  Saffman.^8' 

The  integral  method  may  also  be  applied  to  the  species  conservation 
equations  and  flame  density  equation,  providing  suitable  representation 

for  the  profiles  may  be  constructed.  Those  chosen  were 


<1  = f-|(n)  + C1  g(n) 

(72) 

<2  - f 2 (n ) C2  g(n) 

(73) 

L(n)  = D (f-|  + C1  g)(f2  + C2  g) 

(74) 

where  the  functions  f -j , f2,  and  g are 
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9 = 


(n-j  - n)^  (n  - n2)^ 
(ni  - n2)4 


(77) 
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C1  = C2  = 7 


2 (1  + $)  + 0.81 5 X 


(78) 
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If  we  utilize  the  integral  technique  to  evaluate  the  fuel  consumption 
rate  in  the  mixing  zone,  as  expressed  by  Equation  67,  we  find  that 
approximately 
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valid  for  limited  variations  of  <J>  from  unity.  To  the  same  approximation, 
the  ratio  of  fuel  consumption  to  fuel  inflow  is 


0.164  1 + 


(]  Ul")X<f>  (]  I (1  + $)  + 0.815  a) 


(80) 


33 


5.  TURBULENT  FUEL  JET 


The  fuel  jet  problem  utilizes  the  same  general  equations  for  the  fluid 
dynamics  as  employed  for  the  mixing  zone.  Section  4,  and  because  we  approx- 
imate the  jet  by  a two-dimensional  momentum  source,  the  solution  has  the 
familiar  similarity  property.  The  similarity,  again,  rests  on  the  assump- 
tion that  the  heat  addition  is  small  so  that  the  fluid  mechanical  and 
combustion  problems  separate.  The  combustion  part  of  the  problem  does  not 

have  a similarity  solution,  however,  simply  because  there  is  a finite  amount 
of  fuel  injected  and  this  is  eventually  consumed,  so  that  the  chemical 
reaction  ceases.  Moreover,  the  classical  "point"  jet  similarity  solution 
carries  an  initially  zero  mass  flux  while  the  fuel  jet  problem  requires  a 

fixed  flux  of  fuel  to  be  injected  at  the  origin.  This  difficulty  is 
resolved  in  the  familiar  manner  by  displacing  the  effective  origin  to  an 
appropriate  point  downstream  of  the  momentum  source. 

Because  the  combustion  solution  is  ron-similar,  the  differential 
equations  for  the  chemical  species  and  the  flame  surface  density  do  not 
reduce  to  ordinary  differential  equations,  but  become  partial  differential 
equations  in  terms  of  the  fluid  dynamic  similarity  variables.  Simplifica- 
tion of  the  resulting  problem  is  achieved  by  the  familiar  procedure  of 
constructing  a linear  combination  of  the  various  mass  fractions  in  which 
specific  atom  concentrations  are  conserved  and,  hence,  satisfy  a homogeneous 
differential  equation.  In  the  present  case,  it  proves  possible  to  develop 
such  a combination  of  fuel  and  oxidizer  mass  fractions  that  is  equal  to  the 
similarity  velocity  distribution.  The  mathematical  problem  then  reduces  to 
a pair  of  partial  differential  equations,  one  linear  and  one  non-linear, 
with  an  algebraic  integral  which  allows  determination  of  the  individual 
species  mass  fractions. 

The  similarity  variables  for  the  two-dimensional  jet  follow  from  the 
observations:  i)  that  there  is  no  natural  length  and,  hence,  the  jet  should 

spread  linearly  along  lines  n = x/y  = constant,  and  ii)  that  the  momentum 

o 

injected  by  the  jet  is  conserved  so  that  the  integral  / u dy  over  the  jet 

is  constant.  But  since  the  jet  width  increases  directly  with  x,  the  velocity 

-1  /2 

distribution  must  vary  as  x It  also  follows  that  the  flow  in  the  jet 

1/2 

varies  as  x , so  that  the  injected  mass  flow  is  zero,  the  point  that 
prompted  the  earlier  observation  about  fuel  flow. 
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The  appropriate  variables  are  then 


K = x 
n = y/x 


(81) 


and  if  we  denote  the  momentum  efflux  from  the  jet,  per  unit  depth  of  the 
x-y  plane,  as  py,  where  p is  the  constant  density  of  the  fluid,  then  it  is 
appropriate  to  write  the  stream  function 

>l>  = F(n)  (82) 

e 1/2 

For  the  similarity  to  hold,  it  is  clear  from  Equation  7 that  - = V ^ E,  ' » 

-3/2  03 

and  from  Equation  8 that  to  ^ £ . Therefore,  it  is  appropriate  to  write 

1/2 

= “372'  ^(n)  (83) 
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e = |E(n)  (84) 

with  fi(n)  and  E(n)  dimensionless.  Under  these  definitions,  the  velocity 
components  are 

U = Vf  F'  (85) 

v =Vt[l  F - "F,l  <86) 


and  Equation  6 becomes 


The  momentum  equation  may  be  integrated  directly  to  give 
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where  the  constant  of  integration  vanishes  because  of  the  conditions  on  F' 
and  E/ft  at  the  edge  of  the  jet.  Equation  7 reduces  to  the  form 

-F'  E^FE'=7®{§Ej+“*lF"!E-Ea  (“) 

where,  following  Saffman's  choice,  we  have  taken  8*  = 1 . We  note  further 
by  the  boundary  conditions  and  symmetry  of  the  problem,  | F" | = - F"  for 
n > 0 in  which  region  the  calculation  will  be  carried  out.  Finally, 

Equation  8 becomes 

- 3F'  ft2  - F ft  ft'  = (E  ft1)  + oi|  F"  | ft2  - 8ft3  (89) 


The  conditions  on  the  dependent  variables  are  that,  at  the  edge  of  the 
jet,  F'(n-j)  = ft(n-j)  = E 1 ( n i ) = 0,  while  on  the  symmetry  axis,  F(0)  = E' (0)  = 
ft1  (0)  = 0.  In  addition,  the  momentum  p per  unit  mass  is 


2 


from  which  it  follows  that 


dy  = y 


(90) 


In  treating  the  chemical  composition  and  flame  density,  we  make  cor- 
responding assumptions  that  were  made  for  the  mixing  layer  concerning 
molecular  weights  of  the  constituent  and  the  possibility  of  grouping  an 
inert  diluent  with  any  of  the  other  constituents.  Thus,  we  deal  with  only 
three  species  that  satisfy  equations 
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Now  by  virtue  of  the  linearity  of  these  equations  in  their  respective  con- 
centrations and  the  fact  that  the  consumption  terms  differ  by  only  multip- 
licative constants,  it  is  possible  to  write  linear  combinations  of  the  mass 
fractions  that  are  conserved.  The  obvious  one,  which  we  have  utilized  in 
the  last  example,  is 

K i + *2  + Kg  = 1 (94) 


where  the  normalization  is  Included  in  the  definition  of  mass  fraction.  A 
second  linear  combination,  of  use  to  us  in  the  present  example,  is 


J = K 


1 


+ 


1 + <j>  k3 


(95) 


which  satisfies  the  homogeneous  equation 


and  vanishes  for  large  values  of  y,  outside  the  jet.  Therefore,  because  it 
satisfies  the  same  equations  and  the  same  remote  conditions  as  U(x,y),  it 
is  proportional  to  this  function  and  the  constant  of  porportionality  must 

be  given  by  an  integral  condition,  similar  to  the  momentum  integral. 

Equation  90,  for  the  jet  itself. 

This  condition  can  be  determined  by  stipulating  that  the  jet  initially 
injects  a volume  flow  rate  y-j  of  fuel  which  implies  that  the  entire  injected 
flow  rate  of  fuel  plus  diluent  is  y-j  [1  + k^(o)/k^ (0)].  Far  from  the  point 
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of  injection,  where  the  fuel  is  entirely  reacted  to  combustion  products, 
the  jet  will  carry  a flow  rate  of  combustion  products 


<4(") 

— 1 

Y1 

1 + Kl^_ 

rL 

1 + 


K4H 

k2H 


(97) 


in  which  the  injected  flow  has  been  augmented  by  the  chemically  correct 
amount  of  the  ambient  fluid  to  react  the  fuel  in  the  stoichiometric  ratio. 
Calling 


Y = Y 
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<4(0) 

iojT6T 


(98) 


expression  97  becomes 


1 + 


1 ( k2(°°)  + K^co^\fK-]  (°) 

f yK-j(o)  + k4(o) 


y(1  + <j>) 


(99) 


Thus  the  integral  condition  on  the  reaction  products  is,  where  the  reaction 
is  complete, 

*1 

lim  2 f U k,  dy  = y(l  + <f>)  (100) 


But,  because  ^ 0 for  large  x,  indicating  that  the  fuel  is  eventually 

reacted,  Equation  95  gives 


so  that 


lim  J = -i — r lim  <_ 

X->  00  V X->  00  J 


1 i m 2 f U J dy  ~ y 


(101) 


To  express  this  condition,  it  has  been  necessary  to  introduce  a volume 
flow  rate  y of  material  that  initially  constitutes  the  jet  which,  in  turn, 
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must  carry  the  momentum  y.  This  is  clearly  not  admissible  in  the  theory  of 

the  point  jet,  for  which  y = 0,  because  between  the  flow  rate  and  momentum 

2 

flux,  we  may  define  a length  y /y  and  a velocity  y/y  that  characterize  a 
jet  of  finite  dimensions.  The  appropriate  viewpoint,  however,  is  similar 

to  the  description  of  a finite  jet  by  means  of  a point  jet  with  a displaced 
origin,  and  shares  with  that  approximation  the  difficulty  that  some  early 
portion  of  the  jet  structure  is  described  inaccurately. 

Now  because  the  function  0(x,y)  is  proportional  to  the  velocity  U(x,y), 
it  may  be  written  in  similarity  coordinates 


J 


(102) 


where  d is  an  unknown  length.  Substituting  this  expression  in  Equation  101 
gives 


lim 

£-X» 
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y 


(103) 


But  the  integral  is  Independent  of  S and  holds  for  all  K so  that  the  normal- 
ization of  the  momentum  integral,  Equation  90,  permits  determination  of  the 
unknown  length  scale, 


_ XI 


004) 


which  we  have  defined  earlier  as  the  effective  height  of  the  jet  discharge. 

The  explicit  solution  for  J that  we  have  obtained  permits  us,  using 
Equations  94  and  95,  to  write  the  oxidizer  mass  fraction  in  the  form 

<2  = 1 ■ <i  " (1  + <f>)(J  - <i ) = 1 + 4>  k-|  - (1  + <f>)  J 

and  therefore 


<2  s 1 + <P  <]  - (1  + 4>)  ^rF'(n) 


(105) 
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This  algebraic  integral,  together  with  Equation  94,  makes  it  necessary  to 
determine  only  two  unknowns,  ^(5,0)  and  sU»n)»  by  the  solution  of  differ- 
ential equations. 

To  complete  the  formulation  of  the  jet  problem,  we  shall  consider  the 
expressions  for  fuel  conservation.  Equation  91,  and  the  flame  surface  density, 
Equation  2,  the  latter  written  in  the  form 


U — + V 
8X 


g(+>  foifi 


(106) 


This  form  of  the  flame  density  equation  differs  slightly  from  that  we 

employed  previously,  in  the  form  of  the  flame  shortening  terms.  Here,  the 
two  terms  are  replaced  by  a single  one  that  reproduces  the  general  physical 

idea  and  preserves  the  important  behavior  of  the  term  at  k-|  = 0 and  at 
k2  = 0,  since  the  two  never  vanish  simultaneously. 

It  remains  to  write  Equations  91  and  106  in  similarity  variables  and 
to  make  the  appropriate  choices  for  the  detailed  representations  for  and 
E . Note  that  if  there  were  no  reactant  consumption  terms  in  Equation  91 
for  the  fuel  mass  fraction,  the  quantity  ic^  would  also  be  proportional  to 

the  velocity  and  to  emphasize  this  fact,  we  write 


;1  = Vr  k1U,”) 


(107) 


remembering  that  in  the  absence  of  combustion,  k-|(6,n)  - F'(n).  Further- 
more, it  will  prove  convenient  to  define  a new  variable  (cf.  Equation  59) 
related  to  the  flame  density 


1 


g(<t>) 


E(e,n)  = LU.n) 


(108) 


where  it  is  particularly  to  be  noted  that  a dimensionless  physical  parameter 
has  been  constructed  of  the  molecular  diffusivity  and  the  volume  flow  y per 
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unit  depth.  It  is  then  a matter  of  a straightforward  reduction  to  show  that 
k-]U»ri)  and  LU,o)  satisfies  the  partial  differential  equations 


8k, 


F1  S JL.Ip  k,  -If8"1  - 3 


(,09) 
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and 


pt  r^k_  — pl|_  |pn||_lplk  = l_ 
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(110) 


where  now,  from  Equation  105, 


e2  = 1 


/|[f-  + * (f*  - k^e.n))] 


(in) 


For  use  in  some  of  the  calculations  that  follow,  and  because  of  their 
intrinsic  value,  the  fuel  and  flame  density  equations  may  be  integrated 
across  the  jet  to  give  the  resulting  relations 


"1 


, d 
5 


/ F'  k i d"  = - \f 


/InT  L dn 


(112) 


-n. 


and 


'1  "1  '1 

(53s  -j]f  F'Ldn»«y  |F"jL  dn  - X /[FT  dn  (113) 


'1 


-n 


-n 
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The  single  term  on  the  right-hand  side  of  Equation  112  represents  the 
integrated  fuel  consumption,  while  those  terms  on  the  right-hand  side  of 
Equation  113  represent,  respectively,  the  flame  surface  growth  by  stretch- 
ing and  the  flame  shortening  by  mutual  annihilation. 


The  problem  of  the  turbulent  jet  has  been  treated  numerically  utilizing 
an  integral  technique  and  reasonable  representations  of  the  profiles.  In 
chosing  the  general  profiles  to  represent  the  vorticity  and  energy  densities 

in  the  turbulence,  numerical  calculations  done  by  Dr.  F.  Milinazzo,  privately 
communicated  to  the  authors,  were  very  useful. 
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Again,  omitting  the  tedious  algebraic  details  that  characterize  an 
integral  method,  it  was  found  that  the  velocity,  vorticity,  and  energy 
distributions  could  be  represented  adequately  as 


F'  = F'(0) 


1 + 3 (1  - 6*) 


si  = n(o)  (1  - 6 ) 


E = E(0)  H-.-6.2)  |7  - 3 (1  - 62)2| 


OH) 

(115) 


(116) 


where  F'(0),  ft(0),  and  E (0)  are  the  values  of  these  variables  on  the 
symmetry  axis  and 

« = n/n-j  (117) 


and  their  numerical  values  are 


F'(0) 

fi(0) 

E(0) 

nl 


1.995 

3.570 

0.724 

0.290 


(118) 


As  before,  the  values  of  the  universal  constants  in  the  turbulence  model 
were  taken  as  a = 0.2,  a*  = 0.3,  and  6 = 5/3. 


Proceeding  to  the  combustion  portion  of  the  model,  we  wish  to  reduce 
the  integral  relations,  Equations  112  and  113,  to  ordinary  differential 
equations  by  the  choice  of  fuel  concentration  and  flame  density  profiles. 

The  results  are  made  more  manageable  by  the  separation  of  i and  n depend- 
ence of  these  functions  to  the  degree  reasonable.  With  this  aim  in  mind, 
we  choose  to  represent 

k^f-.n)  » k(0  F'(n)  (119) 


which  for  k(^)  = 1,  yields  an  exact  result  in  the  absence  of  chemical 
reaction,  and,  consequently,  we  know  in  general  that  k(0  < 1.  The  fuel 
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mass  fraction  then  becomes 

*,(«.»)  ■ y|k(e)  F*  (n)  (120) 

so  that,  for  the  physical  bound  that  ^ < 1 to  hold,  there  is  a lower  limit 
on  the  value  of  g for  which  our  representation  is  valid.  Clearly  this  is 


Z 

~3 


min 


3.980 


(121) 


and  we  shall  consider  this  as  the  minimum  value  of  g for  which  we  can 
expect  a reasonable  representation  of  the  jet  flame. 

The  choice  for  the  flame  density  representation  will  be 

L(?,n)  = iti)  k-j (g ,n)  <2(g,n)  (122) 


which  has  the  obvious  advantages  not  only  of  behaving  correctly  at  the  jet 
boundaries  but  in  making  the  second  integral  on  the  right-hand  side  of 
Equation  113  more  convenient  to  handle.  Substituting  the  representations 
for  k-|(g,n)  and  K^U.n),  Equations  111  and  119,  respectively,  we  write 


L(g»n)  = g(0 


1 + <t>  - 4>  k(g) 


(123) 


where  g(g),  formally  equal  to  s,(g)  k(g),  appears  as  the  second  unknown 
function.  It  is  assumed  in  the  representations  of  both  ^(g.n)  and  L(g,n) 
that  F'(n),  where  it  occurs,  is  given  by  Equation  115  and,  hence,  known. 


Again  omitting  much  laborious  detail  and  denoting  /g/d  s z,  the  species 
concentration  integral,  Equation  112,  may  be  written 


z d_  ($k)  = _ 3.0737 


1 - 1.676  ( 

1 + <t>  - 4>k  j 

L v 

2 )\ 

(124) 


where  it  is  convenient  to  use  the  group  <|>k  in  the  numerical  calculations. 
Proceding  in  a similar  manner,  the  integrated  equation  for  the  flame 
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density,  Equation  113,  becomes 


z dz  = M (5M  + 1,2295  N " 2)  + 3-0737  N9 

where  /,  , ±,  \ 

M = 1 - 1.6755  (126) 

N » 1 - 1.4109  (L±±.-.±^  (127) 

The  numerical  integration  requires  specification  of  the  equivalence  ratio 

<j)  and  the  initial  value  k(zQ)  = 1.  In  addition,  a value  for  the  universal 

constant  X must  be  selected  and  the  initial  value  g(z  ) specified.  One  of 

o 

the  difficulties  in  our  point  jet  approximation  is  that  the  initial  develop- 
ment of  the  jet  of  finite  cross  section  is  not  described  correctly.  There- 
fore, the  present  theory  gives  no  description  of  the  transition  between 
the  mixing  zone-dominated  portion  and  the  fully-developed  jet  and,  hence, 

provides  no  information  concerning  the  flame  density  with  which  to  initiate 
the  jet  calculation.  It  is  possible  to  make  an  estimate  for  g (zQ ) from  the 
results  obtained  in  Section  4 for  the  flame  structure  of  the  mixing  zone. 
While  a matching  procedure  has  not  been  carried  out  in  detail  for  various 
values  of  $,  it  appears  that  0.1  < g(z0)  < 0.3  is  reasonable  for  non-extreme 
values  of  the  equivalence  ratio. 

In  the  calculations  that  have  been  carried  out,  g(zQ)  has  generally 
been  taken  as  0.2,  and  the  effects  of  variations  from  this  value  have  been 
explored  only  to  a limited  extent.  The  item  of  principle  interest  in  the 
calculation  is  the  variation  of  the  fuel  concentration  along  the  jet  axis. 

It  is  this  quantity  that  is  frequently  measured  and,  properly  interpreted, 
gives  the  best  information  on  the  burn-out  of  the  jet.  Thus  since  our  fuel 
has  been  approximated  in  the  form 

F'(n)  k(x/d ) 

it  will  suffice  to  show  the  values  of  k(x/d)  and  infer  that  the  centerline 
fuel  mass  fraction, 
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<1  (x,0)  = yf  F * ( 0 ) k(x/d ) 


(128) 


decreases  somewhat  more  rapidly  with  x.  Because  the  fuel  concentration  of 
the  non-reacting  jet  will  decrease  as  /d/x  along  the  axis  [k ( x/d ) =1],  the 
observed  decrease  of  k(x/d)  with  x is  due  only  to  the  chemical  reaction, 
rather  than  to  the  normal  dilution  due  to  mixing. 


Figure  12  shows  the  value  of  k(x/d)  in  terms  of  distance  from  the  fuel 
injection  point  for  three  different  values  of  x.  It  will  be  recalled  that 
X is  the  single  additional  universal  constant  that  enters  in  the  coherent 
flame  model  and  is  associated  with  the  flame  shortening  mechanism.  All 
three  curves  were  computed  for  the  equivalence  ratio,  $ = 1.0.  The  geometry 
of  the  curves  appear  quite  reasonable  when  compared  with  the  experimental 
results  shown  in  Figures  14  through  17  of  Reference  3.  The  effect  of  x is 
quite  clear  from  our  calculations.  For  large  X,  the  flame  shortening 
process  is  accentuated  which,  in  turn,  leads  to  a low  flame  density  and  a 
long  flame.  Smaller  values  of  x produce  a correspondingly  denser  flame 
structure  and  shorter  flames.  As  Figure  12  indicates,  the  results  are 
quite  sensitive  to  X and  here  it  should  be  possible  to  obtain  a reasonably 
good  determination  of  the  value  of  x from  the  fuel  jet  experiments  that  are 
available.  This  systematic  and  detailed  comparison  has  not  yet  been  made, 

but  a superficial  examination  would  indicate  the  X will  be  some  where 
between  0. 1 and  0.5. 


Figure  13  is  presented  to  show  the  effect  of  equivalence  ratio,  4>, 
on  the  flame  length  when  the  value  of  x has  been  fixed.  Clearly  large 
equivalence  ratios  lead  to  long  flames  and  low  equivalence  ratios  lead 
to  shorter  ones.  This  behavior  results  from  the  relatively  larger  amount 
of  oxidizer  that  must  generally  be  entrained  and  mixed  for  a flame  with 
high  equivalence  ratios  than  for  one  of  low  equivalence  ratio.  Familiar 
examples  of  high  and  low  equivalence  ratio  are  the  methane-air  flame  and 
the  hydrogen-air  flame,  respectively.  The  results  of  Figure  13,  computed 
using  the  value  x = 0.5,  suggest,  upon  comparison  with  the  experiments, 
that  a somewhat  smaller  value  of  x would  be  appropriate. 
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Figure  12.  Effect  of  Flame  Shortening  Constant  on  Length  of  Turbulent  Diffusion 
Flame  from  Fuel  Jet.  6 = 1.0. 


Effect  of  Equivalence  Ratio  on  Length  of  Turbulent  Diffusion 
Flame  from  Fuel  Jet.  X = 0.5. 
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JO-  ABSTRA^  (Continue  on  reveree  side  If  neceeeery  end  Identity  by  block  number) 

A description  of  the  turbulent  diffusion  flame  is  proposed  in  which 
the  flame  structure  is  composed  of  a distribution  of  laminar  diffusion 
flame  elements,  whose  thickness  is  small  in  comparison  with  the  large 
eddies.  These  elements  retain  their  identity  during  the  flame  development; 
they  are  strained  in  their  own  plane  by  the  gas  motion,  a process  that  not 
only  extends  their  surface  area,  but  also  establishes  the  rate  at  which  a 
flame  element  consumes  the  reactants.  Where  this  flame  stretching  process 
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has  produced  a high  flame  surface  density,  the  flame  area  per  unit  volume, 
adjacent  flame  elements  may  consume  the  intervening  reactant,  thereby 
annihilating  both  flame  segments.  This  is  the  flame  shortening  mechanism 
which,  in  balance  with  the  flame  stretching  process,  establishes  the  local 
level  of  the  flame  density.  The  consumption  rate  of  reactant  is  then  given 
simply  by  the  product  of  the  local  flame  density  and  the  reactant  consump- 
tion rate  per  unit  area  of  flame  surface.  The  proposed  description  permits 
a rather  complete  separation  of  the  turbulent  flow  structure,  on  one  hand, 
and  the  flame  structure,  on  the  other,  and  in  this  manner  permits  the 
treatment  of  reactions  with  complex  chemistry  with  a minimum  of  added 
labor.  The  structure  of  the  strained  laminar  diffusion  flame  may  be 
determined  by  analysis,  numerical  computation,  and  by  experiment  without 
significant  change  to  the  model 

The  flame  density  and  the  mass  fractions  of  reactant  are  described 
by  non-linear  diffusion  equations  in  which  those  equations  for  the  reactants 
each  contain  a consumption  or  production  term  proportional  to  the  local 
flame  density.  The  flame  density  equation  contains  a production  term 
associated  with  flame  surface  stretching  and  a consumption  term  describing 
the  flame  shortening  by  mutual  annihilation.  Each  of  the  equations  contains 
a turbulent  diffusion  term  utilizing  a scalar  diffusivity.  The  model  of 
inhomogeneous  turbulence,  proposed  by  Saffman,  completes  the  description^- 
of  the  problem  and  couples  with  the  flame  and  composition  equations  to 
determine  the  velocity  distribution  and  the  turbulent  diffusivity.  A 
single  additional  universal  constant,  over  those  appearing  in  Saffman's 
model,  is  required  in  the  model  equations  for  the  flame. 

The  coherent  flame  model  has  been  applied  to  diffusion  flame  structure 
in  the  mixing  region  between  two  streams  and  predicts  correctly  the  result 
that  the  reactant  consumption  per  unit  length  of  flame  Is  independent  of 
the  distance  from  the  initiation  of  mixing.  In  this  example  which  Is 
carried  out  for  small  density  changes,  both  the  fluid  mechanical  and  flame 
variables  possess  similarity  solutions. 

The  coherent  flame  model  is  also  applied  to  the  turbulent  fuel  jet 
which  clearly  does  not  have  a similarity  solution  simply  because  the  finite 
mass  flow  of  fuel  is  eventually  consumed.  The  problem  is  solved  utilizing 
an  integral  technique  and  numerical  integration  of  the  resulting  differential 
equations.  The  model  predicts  the  flame  length  and  superficial  comparison 
with  experiments  suggest  a value  for  the  single  universal  constant.  The 
theory  correctly  predicts  the  change  of  flame  length  with  changes  In  stoich- 
iometric ratio  for  the  chemical  reaction. 
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